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ABSTRACT 

We present an early broad-brush analysis of HerschelfPACS observations of star-forming galaxies in 8 galaxy clusters drawn from our survey of 30 
clusters at z=:0.2. We define a complete sample of 192 spectroscopically confirmed cluster members down to Ltir>3x lO'^L© and Lk>0.25Lq. The 
average A'-band and bolometric infrared luminosities of these galaxies both fade by a factor of ~2 from clustercentric radii of ~2r2oo to ~0.5r2oo, 
indicating that as galaxies enter the clusters ongoing star-formation stops first in the most massive galaxies, and that the specific star-formation 
rate (SSFR) is conserved. On smaller scales the average SSFR jumps by ~25%, suggesting that in cluster cores processes including ram pressure 
stripping may trigger a final episode of star-formation that presumably exhausts the remaining gas. This picture is consistent with our comparison 
of the //er.st77e/-detected cluster members with the cluster mass distributions, as measured in our previous weak-lensing study of these clusters. 
For example, the spatial distribution of the Herschel sources is positively correlated with the structures in the weak-lensing mass maps at ~5cr 
significance, with the strongest signal seen at intermediate group-like densities. The strong dependence of the total cluster IR luminosity on cluster 
mass - LtirocM^j^j^i - is also consistent with accretion of galaxies and groups of galaxies (i.e. the substructure mass function) driving the cluster IR 
luminosity. The most surprising result is that roughly half of the Herschel-AeXe.cXe,A- cluster members have redder 5 ioo/'S24 flux ratios than expected, 
based on the Rieke et al. models. On average cluster members are redder than non-members, and the fraction of red galaxies increases towards the 
cluster centers, both of which indicate that these colors are not attributable to systematic photometric errors. Our future goals include to intepret 
physically these red galaxies, and to exploit this unique large sample of clusters with unprecedented multi-wavelength observations to measure the 
cluster-cluster scatter in SO progenitor populations, and to intepret that scatter in the context of the hierarchical assembly of clusters. 
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1. Introduction 

Lenticular galaxies (hereafter SOs) are mainly found in the cores 
of galaxy clusters at low redshift (e.g. Dressier et al. 1997; Smith 
et al. 2005; Postman et al. 2005). There is a broad consensus that 
they are the descendants of gas rich spiral galaxies that have been 
accreted from the surrounding filamentary structure. However 
the physics of how spirals are transformed into SOs remains 
largely unconstrained, with numerous "SO progenitor" popula- 
tions (e.g. Moran et al. 2006; Poggianti et al. 2000; Geach et al. 
2006; Haines et al. 2009a - hereafter H09a) and physical pro- 
cesses (e.g. Gunn & Gott 1972; Moore et al. 1999) discussed in 
the literature. 

The broad range of cluster-centric radii at which various SO 
progenitors are found reflects the fact that different physical pro- 
cesses act in different environments, for example ram pressure 
stripping is more effective closer to cluster centers where the 
intracluster medium (ICM) is denser, and galaxy-galaxy merg- 
ing is more effective in galaxy groups that are falling into the 
cluster than in the cluster cores. Moreover, the observational 
signatures of SO progenitors are diverse, ranging in wavelength 
from ultraviolet (UV) emission from A stars in galaxies whose 
star-formation (SF) has been recently quenched, through optical 
spectral features including Balmer absorption lines, to mid/far- 
infrared (IR) emission from dust heated by SF (e.g. Moran et al. 
2007; Poggianti et al. 2000; Haines et al. 2009b). 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



Mid- and far-IR properties of cluster galaxies have been 
studied previously with IRAS (e.g. Leggett et al. 1987; Doyon 
& Joseph 1989), ISO (see Metcalfe et al. 2005 for a review), 
and Spitzer (e.g. Geach et al. 2006; Fadda et al. 2008; Haines et 
al. 2009a,b; Bai et al. 2009). A key result from these IR stud- 
ies is that a significant fraction of the total SF in galaxy clusters 
is obscured by dust. The inferred levels of SF naturally fit the 
hypothesis that bulge dominated SOs are descended from late- 
type spirals. It has also been suggested that dusty SO progenitors 
are more common in dynamically active, i.e. merging, galaxy 
clusters than in so-called "relaxed" clusters (e.g. Metcalfe et al. 
2005; Geach et al. 2006; Miller et al. 2006). However it has thus 
far been difficult to test this idea robustly because the intrinsic 
scatter in levels of SF in clusters appears to be large, (as noted 
by Kodama et al. 2004), and the sample sizes observed to date 
are small (i.e. < 2) within any given redshift bin - although see 
H09a for a recent counter-example. 

We are therefore conducting a systematic wide-field survey 
of a large statistically well-defined sample of galaxy clusters in 
a narrow redshift slice at z-0.2, as part of the Local Cluster 
Substructure Survey (LoCuS^. Our goals are to compile a 
complete inventory of SO progenitors using data from the far- 
UV to far-IR, and to relate these populations to the underlying 
gas physics and hierarchical structure of the host galaxy clusters. 
We aim to delineate the different physical processes responsible 
for galaxy transformation in clusters and their suiTounding large 
scale structure, and thus constrain the amplitude of the different 
physical pathways from spiral to SO morphology, and how these 
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relate to the dynamical state of the clusters. Our Open Time Key 
Programme observations with Herschel (Pilbratt et al., 2010), 
supplemented by existing Spitzer mid-IR observations provide 
the all-important measurements of the bolometric IR luminosity 
and mid/far-IR colors of dust-reddened/obscured SO progenitors. 
We assume //o=70km s"', Qm=0.3, Qa=0.7. In this cosmol- 
ogy Ikpc at z=0.2, subtends 0.3". All cluster masses and radii 
relative to an over-density are derived from the weak-lensing 
analysis of Okabe et al. (2010; hereafter OklO). 



2. Survey Design 

Our survey goals include understanding the physical reasons for 
the large cluster-cluster variations in SF rate (SFR), and the full 
range of physical processes responsible for transforming spiral 
galaxies into SOs. We therefore require a large sample of clusters 
in order to sample thoroughly the underlying cluster population 
and the various SO progenitor populations that they host. Our 
sample of 30 clusters therefore will allow us to study, for exam- 
ple, how total integrated cluster SFRs depend on global cluster 
properties such as cool core strength, and substructure fraction, 
in ~3-6 cluster bins with ~5-10 clusters per bin. Based on pre- 
vious IR and UV studies of SO progenitors, we expect ~30-50 
such objects per cluster Our sample of 30 clusters should there- 
fore deliver a sample of >1000 SO progenitors. 

With current observing facilities and observed cluster sam- 
ples, this study is prohibitvely expensive at high redshift because 
of the requirement for wide-field and moderately deep data on a 
large sample. We therefore concentrate on clusters at z^O.2; at 
this redshift gravitational lensing is an efficient probe of the dark 
matter distribution in the cluster in-fall regions using Suprime- 
CAM on the Subaru telescope, and yet follow-up Herschel ob- 
servations of a sample of 30 clusters are feasible. 

The cluster sample is a subset of those in the ROSAT All- 
sky Survey catalogs (Ebeling et al. 1998, 2000; Bohringer 
et al. 2004) that satisfy the following criteria: 0.15<z<0.3, 
«H<7xlO^"cm"^, and -70°<(5<+70°, and that were observable 
with Subaru/Suprime-CAM on the nights assigned to us (OklO). 
The sample is therefore blind to the thermodynamic, and hier- 
archical assembly history of the clusters, other than the use of 
X-ray luminosity as a proxy for mass-selection. The distribu- 
tion of the X-ray luminosities of clusters in our sample is sta- 
tistically indistinguishable from that of a volume-limited sam- 
ple satisfying Lx £'(z)"^'^>4.2xl0'*'*erg s"', where the scaling of 
Lx with E{z)^-Q.M {^+z)^+0.a approximates mass-selection fol- 
lowing Popesso et al. (2005) - see OklO for more details. 

Our multi-wavelength observations of this sample (§3) span 
at least a clustercentric radius of 12.5arcmin on the sky, equating 
to ~ 1 .5r2oo for a typical cluster in our sample. This physical field 
of view is sufficiently large to probe all of the physical processes 
expected to play a role in galaxy transformation (Fig. 1). 



3. Observations and Data Analysis 

The 8 clusters discussed in this letter were observed with 
the Photodetector Array Camera and Spectrometer (PACS; 
Poglitsch et al. 2010) across a 25'x25' field of view at 100 
and 160yum in scan map mode at 20arcsec/s in November and 
December 2009. Each cross scan was repeated 7 times, giving a 
total exposure time of 93.2s/pix. We first processed the data us- 
ing standard hipe routines (Ott 2010). Then all sources detected 
at >2.5cr in the first pass reduced data were masked using 15" 
circular apertures, and the data were high pass filtered with a 
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Fig. 1. This figure shows simple physically well-motivated mod- 
els for the relative strength of ram pressure stripping, harass- 
ment, and galaxy-galaxy merging in the eight clusters discussed 
in this letter. The vertical black dotted lines show the stripping 
radius for each cluster - the radius within which ram pressure 
is able to strip a disk galaxy with rotation speed of 160 kms"'. 
This calculation combines Sanderson et al.'s (2009) Chandra- 
based gas density profiles with OklO's weak lensing-constrained 
NFW density profiles, and assumes that the typical radial veloc- 
ity of infalling galaxies is 0.9yc at the virial radius where Vc 
is the circular velocity of the cluster dark matter halo, follow- 
ing the analytic model of McCarthy et al. (2008). Red dashed 
curves show the relative efficiency of harassment, based on the 
harassment rate of /h oc pgair^ (Moore et al. 1998), assuming 
that the total mass distribution traces the galaxy density (Lin 
et al. 2004) - i.e. we substitute pnfw from OklO for pgai. The 
galaxy-galaxy merger rate is traced by the solid blue curves, 
which show the radial number density profile of galaxies that in- 
habit group-like environments, i.e. a 3-dimensional galaxy den- 
sity of 2<p3D<10Mpc"^, derived from our near-IR photometry 
(H09a) and large spectroscopic catalogs (Hardegree-UUman et 
al. in prep.). 



filter 25 and 30 frames wide at lOO/^m and 160;um respectively. 
The final maps were then constructed using the photproject rou- 
tine. The angular resolution of the final reduced frames is 6.8" 
at lOOyum and 1 1 .4" at 160/im. 

Sources were extracted using SExtractor, employing circu- 
lar apertures of 12" and 16" diameter at 100 and 160jum. The 
point spread function in the 100/zm frame is slightly de-graded 
from that derived from the PACs calibration maps of Vesta, due 
to imperfect spatial calibration within scans where bright sources 
cannot be used to register the individual exposures. We therefore 
applied an empirical aperture correction of 1.65+0.03, based 
on brightest isolated sources in the 100/im maps. The standard 
Vesta aperture correction of 1.675 was applied at 160;um. The 
90% completeness limits are 18mJy at lOOyum and 28mJy at 
160yum. Calculation of the total IR luminosities (8</l<1000/im) 
discussed in §4 are described by Haines et al. (2010). 

We also use our wide-field data from 5/?/fzer/MIPS 
(24/zm; H09a), UKIRT/WFCAM (///T-bands; H09a), Chandra 
(Sanderson et al., 2009), and Subaru/Suprime-CAM (V//' -bands; 
Okabe & Umetsu 2008; OklO). We have spectroscopically iden- 
tified 92% of the sources down to the 100 /im detection thresh- 
old using MMT/Hectospec (Fabricant et al. 2005), in addi- 
tion to securing ~200-300 cluster galaxy redshifts per cluster 
(Hardegree-Ullman et al., in prep.). 
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Fig. 2. A'-band luminosity (Lk) versus total infrared luminos- 
ity (Ltir) for the Herschel-detected cluster members, plotted as 
open circles. The horizontal and vertical solid lines show the 
completeness limits of our MMT/Hectospec spectroscopic sur- 
vey and our Herschel photometry respectively; the horizontal 
and vertical dot-dashed lines show L^ from Lin et al. (2004) and 
the definition of LIRGs respectively. The diagonal dashed line 
shows LkocLtir, and guides the eye for the radial dependence of 
{Lk) and (Ltir) respectively, as shown by the color-coded filled 
circles - see legend. 



4. Results 

We identify 192 Herschel sources down to Ltir>3x1O"^L0 
and Lk>Q-25L^ (Fig. 2) within a clustercentric radius of 
R<l.5r2Qo and lying inside the caustics in the velocity -radius 
plane (see Haines et al. 2010 for an example). Just (24+3)% 
(46/192) of these galaxies are LIRGs (Ltir>1O"L0) and 
none are ultra-luminous IR galaxies (ULIRGs; Ltir>1O'^L0). 
The most luminous galaxy is the brightest cluster galaxy 
(BCG) in A 1835; we also detect the BCG in A2390 (Edge 
et al. 1999; Egami et al. 2006). The typical galaxy has 
La:-1-5x1O"L0 and Ltir-6x1O'°L0. On average, the most 
IR-luminous galaxies are found at projected cluster-centric 
radii of ~1.5r2oo; at larger and smaller radii the average IR- 
luminosity declines. This trend is mirrored by a decline in 
the average A'-band luminosity such that (Ljs^/Lk) is con- 
served down to ~0.5r2oo, interior to which the average IR- 
luminous galaxy jumps from (LTiR/iA:)(^>0.5r2oo)-0.39+0.03 
to (LxiR/ijf >(-R<0.5r2oo)-0.53+0.03 (Fig. 2). 

The 5 100/5' 24 colors of the IR-detected galaxies have a 
prominent excess of flux at lOOyKm relative to that expected 
from the commonly used Rieke et al. (2009) SED templates 
(Fig. 3). Adopting S ioo/524>25 as defining this unexpected red 
population, we find that (63+4)% (121/192) of the spectroscop- 
ically confirmed members are "red", with (78;^g)% (36/46) of 

LIRG members, and (58!^)% (85/146) of sub-LIRG members 
being "red" respectively. The predominance of red LIRGs is 
partly expected given the predicted relationship between lumi- 
nosity and color, however the observed LIRGs lie almost ex- 
clusively red-ward of the Rieke et al. models. The fraction of 
galaxies with red colors also shows a gentle increase towards 



the cluster centers: f^ed'^R^^^'^'^'^'^^ ■ The mean color of Herschel- 
detected cluster members, (S ioo/'S24)=29. 2+0.1, is also slightly 
redder than the mean color of Herschel-detected non-members 
(defined as lying outside the caustics, but within 0.15<z<0.3), 
{S 100/5 24)=27.0±0.1. Both of these correlations point to a phys- 
ical origin for these red colors, rather than systematic photomet- 
ric errors. Indeed, Rawle et al. (2010) find a similar population 
in the core of the Bullet cluster However the small difference 
between the color distributions of members and non-members 
suggest that cluster physics may not be the main determinant of 
the IR colors (see also Pereira et al. 2010). 

The IR luminous population traces the structure of the mass 
distribution of the clusters and the surrounding filamentary struc- 
ture, as traced by OklO's weak-lensing mass maps (Fig. 4). This 
is quantified using the number density of galaxies (normalized 
to the mean galaxy number density in each case) as a function 
of projected cluster mass density obtained from the mass maps 
(Fig. 4). A spatially random distribution of galaxies is consis- 
tent with unity; positive correlation between galaxies and mass 
is greater than unity, and anti-correlation is less than unity. We 
detect correlation between the spatial distribution of the IR lumi- 
nous cluster members and the cluster mass distributions at ~5cr 
significance. 

We also combine the Herschel data with OklO's weak- 
lensing analysis to construct the first ever mass-LxiR relation 
for galaxy clusters (Fig. 3), obtaining Ljs^(<ryi^^\)ccM"^^,^^^, with 
Q'=1.9;^Qg. The fit was done taking into account errors in both 
variables, and was repeated 10"* times, each time drawing 8 clus- 
ters at random with replacement; the error quoted on the slope is 
dominated by the scatter between these bootstrap samples. The 
mass-to-light ratio of a ~1O'^M0 cluster is ~1O^M0/L0, how- 
ever the scaling of mass with far-IR luminosity is inconsistent 
with a constant mass-to-light ratio at ~2cr. 

Cluster IR luminosity density profiles, normalized to the un- 
derlying dark matter density profile from Ok 10, are generally 
increasing functions of clustercentric radius out to at least r2oo 
(Fig. 3). IR luminosity in excess of that expected from a flat 
mass-to-light ratio profile is therefore found at large cluster- 
centric radii, particularly at /?~l-1.5rviiiai, where some of the 
profiles show a pronounced peak. We also note that strong- 
lensing clusters (A 1835, A 1689, A 2219, A 2390) tend to have 
steeper mass-normalized luminosity density profiles than non- 
strong-lensing clusters, three of which are well-known merging 
clusters (A 1763, A 1758, A 1914) in which the merger axis is 
likely close to the plane of the sky. Cluster geometry, e.g. pro- 
late shape and/or merger aligned with the line of sight (strong- 
lensing clusters) versus aligned in the plane of the sky, therefore 
may complicate the detailed interpretation of the luminosity den- 
sity profile shapes. 



5. Summary and Discussion 

We have presented an initial broad-brush analysis of 
HerschelfPACS observations of 25% of our sample of 30 
galaxy clusters at z-0.2, and combined these data with our 
existing Spitzer, Subaru, Chandra, UKIRT, and MMT data. 
The main analysis concentrates on a sample of 192 spectro- 
scopically confirmed cluster members with Ltir>3x1O'"L0, 
Lk>0-25L^, R<l.5r2Qo. The average /T-band luminosity of 
these galaxies fades by a factor of almost 2 from the cluster 
outskirts (~l-2r2oo) to the cluster cores (<0.5r2oo), although the 
average specific star- formation rate, as probed by (Ljiji,) / (Lk) , 
is constant across most of this radial range (~0.5-2r2oo), before 
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Fig. 3. Left - Far-IR color-luminosity relation; the color-luminosity relation predicted by the Rieke et al. (2009) SED templates is 
shown as the black solid line. The vertical and horizontal dashed lines mark the definition of LIRGs, and the nominal color cut 
discussed in §4, respectively. Cluster members are shown as filled black circles, and non-members as open purple circles. Center 
- Cluster virial mass versus bolometric infrared luminosity within the virial radius; the best-fit relation is shown as a solid line. 
Right - Total IR luminosity density profile relative to total mass density profile (§4). Strong-lensing clusters (Richard et al. 2010) 
are plotted as solid curves and non-strong-lensing clusters as dashed curves. The legend lists clusters from most (A 1835) to least 
(Z7160) massive. 
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Fig. 4. Left - Weak lensing mass maps (contours, spaced at Icr mass reconstruction error; OklO; Okabe & Umetsu 2008) in order 
of decreasing cluster virial mass - top left-to-right: A 1835, A 1763, A 1689, A2219; bottom: A2390, A 1914, A 1758, Z7160. 
Spectroscopically confirmed //eric/ieZ-detected members (selected as described in §4) are marked as filled red circles that scale 
logarithmically with Ljir/Lk- North is up and East is left, and the axis tick-marks are spaced at 1 arcmin intervals. Right - Number 
density of galaxies relative to the mean number density of each population as a function of total projected cluster mass density. 
Positive correlation of the i/e^c/ieZ-detected cluster members with the mass distribution in/surrounding the clusters is detected at 
~5cr. The absence of red triangular data points at 2niass>1000MQkpc"^ is because we do not find any non-BCG IR luminous cluster 



members at such high densities. The dip in the signal for cluster members at ~125OM0kpc 
spectroscopic catalog in high density environments due to MMT/Hectospec fiber collisions. 



is due to the incompleteness of our 



jumping by 25% on smaller scales. This suggests that as gas 
rich galaxies fall into the clusters (typically in groups - Fig. 1) 
ongoing star-formation stops first in the most massive galaxies. 
As galaxies reach the cluster cores physical processes that 
operate in high density environments, for example ram pressure 
stripping and harrassment, then appear to trigger a final episode 
of star-formation that presumably exhausts the remaining gas 
supply. 

This picture is consistent with our comparison of the 
//eric/ieZ-detected cluster members with the cluster mass distri- 
butions, as probed by Okabe et al.'s (2010) weak-lensing analy- 
sis. First, the spatial distribution of the Herschel sources is pos- 
itively correlated with the structures in the weak-lensing mass 
maps at ~5cr significance, with the strongest signal seen at inter- 
mediate, group-like densities. Second, the strong dependence of 
the total cluster IR luminosity on cluster mass (Ltir°c Af virial-* ^^ 



consistent with accretion of galaxies and groups of galaxies driv- 
ing the cluster IR luminosity. This is because, assuming that IR 
galaxy mass-to-light ratios are independent of the cluster mass, 
the scaling relation can be understood as stemming from the 
M^ dependence of the substructure mass function seen in the- 
oretical models (e.g. Taylor & Babul 2005). Third, the IR lumi- 
nosity density profiles of the clusters generally increase to large 
radii, with some clusters showing a peak at ~r2oo- This is qual- 
itatively consistent with the enhanced star-formation rates seen 
in in-falling galaxy populations by Moran et al. (2005). 

The most surprising result is that roughly half of the 
Herschel-detected cluster galaxies have excess flux at 100/vm 
over that predicted from current SED models. Cluster members 
are redder than non-members, and we find a shallow trend of in- 
creasing fraction of red IR galaxies towards the cluster centers, 
both of which suggest that this is a physical effect and not caused 
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by systematic photometric uncertainties. This result will be the 
focus of more detailed future investigation. 

Finally, we note that, contrary to previous speculation in the 
literature, we do not find a strong relationship between cluster 
IR luminosity and cluster dynamical state. Observations of the 
full sample will allow us to investigate this issue in more detail. 
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